Vortex Origin of Tricritical Point in Ginzburg-Landau Theory 
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Motivated by recent experimental progress in the criti- 
cal regime of high-T c superconductors we show how the tri- 
critical point in a superconductor can be derived from the 
Ginzburg-Landau theory as a consequence of vortex fluctu- 
ations. Our derivation explains why usual renormalization 
group arguments always produce a first-order transition, in 
contrast to experimental evidence and Monte Carlo simula- 
tions. 

1. The critical regime of old-fashioned superconductors 
can be describes extremely well by the Ginzburg-Landau 
theory [1] at the mean-field level [2,3]. The reason is 
the smallness of the Ginzburg temperature interval ATq 
around the mean-field critical temperature T^ F where 
fluctuation become important [4]. This made discus- 
sions of the order of the superconductive phase transition 
started by Halperin, Lubensky, and Ma in 1972 [5] rather 
academic. 

The situation has changed with the advent of mod- 
ern high-T c superconductors temperature superconduc- 
tors. In these the Ginzburg temperature interval is large 
enough to study field fluctuations and critical behavior. 
Several experiments have found a critical point of the 
XY universality class [6] . In addition, there seems to be 
recent evidence for an additional critical behavior asso- 
ciated with the so-called charged fixed point [7] . In view 
of future experiments, it is important to understand pre- 
cisely the nature of critical fluctuations. 

The Ginzburg-Landau theory [1] describes a supercon- 
ductor with the help of an energy density 

Hty, W, A, VA) = i{[(V - ieA)^} 2 + T\^\ 2 + ||V| 4 } 

(1) 
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where tp{x), A.{x) are pair field and vector potential, re- 
spectively, and e is the charge of the Cooper pairs. The 
parameter r = T/T C MF — 1 is a reduced temperature 
measuring the distance from the characteristic temper- 
ature T^ F at which the r-term changes sign. The the- 
ory needs gauge fixing, which is usually done by setting 
ip(x) = p(x)e l0lyX \ rewriting the covariant derivative oft/' 
as 



Dtp= [i(V0-eA)p + Vp]e 



and eliminating the phase 
gauge transformation A — s 
H(il>, VV>, A, VA 



variable 6{x) 
A + V6/e. 
to the form 



(2) 
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where V(p) is the potential of the p-field: 

V(p) = -/ + 9 -/. 



(3) 



(4) 



The last term in (3) is the famous Meissner-Higgs mass 
rriA = pe [2,3] of the vector potential A, whose analog in 
the gauge theory of electroweak interactions to explain 
why interactions are so much weaker than electromag- 
netic interactions. 

At the mean-field level, the energy density (3) de- 
scribes a second order phase transition. It takes place if 
r drops below zero where the pair field ip(x) acquires the 
nonzero expectation value {ip{x)) = pa = x/—i~/g, the or- 
der parameter of the system. The p-fluctuations around 
this value have a coherence length £ = 2t. The 

Meissner-Higgs mass term in (3) gives rise to a finite pen- 
etration depth of the magnetic field A = \Jvha = 1/poe. 
The ratio of the two length scales k = A/\/2£, which for 
historic reasons carries a factor y2j is the G inzbu rg pa- 
rameter whose mean field value is kmf = yV?. Type 
I superconductors have small values of k, type II super- 
conductors have large values. At the mean-field level, the 
dividing line is lies at k = 1/ V2. 

In high-T c superconductors, field fluctuations become 
important. These can be taken into account by calculat- 
ing the partition function and field correlation functions 
from the functional integral 



Z = / VppVAe 



J d 3 xH% 



(5) 



(in natural units with fc^T =1). So far, all approx- 
imations to Z found since the initial work [5] have had 
notorious difficulties in accounting for the order of the su- 
perconductive phase transition. In [5], simple renormal- 
ization group arguments [8] in 4 — e dimensions suggested 
that the transition should be of first order. The techni- 
cal signal for this is the nonexistence of an infrared-stable 
fixed point in the renormalization group flow of the cou- 
pling constants e and g as a function of the renormaliza- 
tion scale. Due to the smallness of the Ginzburg interval 
ATq, the first order was never verified experimentally. 
Since then, there has been much work [9] trying to find 
an infrared-stable fixed point in higher loop orders or 
by different resummations of the divergent perturbations 
expansions. 

The simplest argument suggesting a first-order nature 
of the transition arises at the mean-field level of the pair 
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field p as follows: The fluctuations of the vector potential 
are Gaussian and can be integrated out in (5). Assuming 
p to be smooth, this may be done in the Thomas-Fermi 
approximation [10], leading to an additional cubic term 
in the potential (4), which becomes 



V(p) 



if + Y 



2tt' 



(6) 



The cubic term generates, for t < c 2 /4g, a second mini- 
mum at 



If r decreases below 



ri = 2c 2 /9g. 



(7) 



(8) 



the new minimum lies lower than the one at the origin 
(see Fig. 1), so that the order parameter jumps from zero 
to 

Pi = 2c/3g (9) 
in a phase transition. At this point, the coherence length 



of the /^-fluctuations £ = 1/yr + 3gp 2 — 2cp has the fi- 
nite value (the same as the fluctuations around p = 0) 



The phase transition is therefore of first-order. 



(10) 



V(p) 
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FIG. 1. Potential for the order parameter p with cu- 
bic term. A new minimum develops around p\ causing a 
first-order transition for t — t±. 

This conclusion is reliable only if the jump of po is suffi- 
ciently large. For small jumps, the mean-field discussion 
of the energy density (6) cannot be trusted. The place 
where the transition becomes second order has, so far, 
never been explained satisfactorily within the Ginzburg- 
Landau theory. An explanation has been found only 
with the help of a dual disorder field theory derived from 
the Ginzburg-Landau theory in Ref. [11,12]. This the- 
ory is constructed in such a way that its Feynman dia- 
grams are direct pictures of the vortex lines of the su- 
perconductor. The dual disorder field theory shows that 



there is a first-order transition only if the Ginzburg pa- 
rameter k = A/\/2£ is smaller than the tricritical value 
K tnc 0.8/ v^. This point is close to the mean-field value 
k = l/y/2 where the superconductor changes from type 
II to type I, and the average short-range repulsion be- 
tween vortex lines changes into an attraction, so that the 
quartic term in the dual field theory becomes negative 
[11]. 

In contrast to the Ginzburg-Landau theory, the vector 
potential of the disorder field theory is massive from the 
outset, so that its fluctuations do not generate a cubic 
term. Instead, they change the sign of the quartic term, 
making the transition first-order for k < K trlc , while leav- 
ing it second-order for k > K trlc . 

The purpose of this note is to show how the tricritical 
point can be derived from the original Ginzburg-Landau 
theory by a proper inclusion of fluctuation corrections. 
The mistake in the above argument lies in the neglect 
of vortex fluctuations. In fact, the transformation of the 
covariant derivative Dtp to the p-6 expression in Eq. (2) 
is false. Since 0(x) and #(x) + 27r are physically indistin- 
guishable — the complex field ip(x) is the same for both 
— the correct substitution is 



£>V = [*(V0 - 2ir0 v - eA)p + V p]e l6 



(11) 



The cyclic nature of the scalar field 9(x) requires the 
presence of a vector field 6 V (x) called vortex gauge field. 
This field is a sum of 5-functions on Voltcrra surfaces 
across which 0(x) has jumps by 27r. The boundary lines 
of the surfaces are vortex lines. They are found from the 
vortex gauge field 6 V (x) by forming the curl 



vx^(x)=r(x), 



(12) 



where j^x) is the vortex density, a sum over 5- functions 
along the vortex lines <5(L;x) = f L dx<5(x — x). Vortex 
gauge transformations are deformations of the surfaces at 
fixed boundary lines which add to v (x.) pure gradients 
of the form V<5(V;x), where 5(V;x) = J v d 3 xS(x — x) 
are (^-functions on the volumes V over which the surfaces 
have swept. The theory of these fields has been devel- 
oped in the textbook [12] and the Cambridge lectures 
[13]. Being a gauge field, 0"(x) may be modified by a 
further gradient of a smooth function to make it purely 
transverse, V • #t( x ) — 0, as indicated by the subscript 
T. Since the vortex gauge field is not a gradient, it cannot 
be absorbed into the vector potential by a gauge trans- 
formation. Hence it survives in the last term in Eq. (3), 
and the correct partition function is 



V0 V T / XVpPAexp 



1 



(V x A) 2 



p 2 e 2 



{A-2TiO v T /ef 
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(13) 



The symbol J T>Q V T does not denote an ordinary func- 
tional integral. It is defined as a sum over all numbers 
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and shapes of Volterra surfaces S in Q V T , across which the 
phase jumps by 2tt [13]. 

The important observation is now that the partial par- 
tition function contained in (13) 



Zi[p] = J V6 v T VAexp 



\ J d 3 x(V x A) 2 

J d 3 x[eA-27r6>^] 2 | (14) 



can give rise to a second-order transition of the XY- 
model type if the Ginzburg parameter k is sufficiently 
large. To see this we integrate out the A-field and obtain 



Z 1 [p] = cxp 
~4w 2 p 2 



p 3 3 

d 3 x 6 9 



6tt 



J ve v T 



(15) 



x exp 
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d 3 x 



1 



-e 



,v 1 



0% 



p 2 e 2 



2 J \2 1 ~+ _ v -' + p 2 e 2- 
The second integral can be simplified to 
4ttV f „ /„„ -V 2 



:0\ 



d 3 a 



-V 2 + p 2 e 2 



.01 



(16) 



Integrating this by parts, and replacing Vi#y Vi<?^ by 
(V x V T ) 2 = j v2 , the partition function (15) without the 
prefactor takes the form 



[p]=JvO v T exp -^f- J d 3 x (? 



1 



-V' z + p 2 e 2 " 



(17) 



This is the partition function of the grand-canonical en- 
semble of closed fluctuating vortex lines. The interaction 
between them has a finite range equal to the penetration 
depth A = 1/pe. 

It is well-known how to compute pair and magnetic 
fields of the Ginzburg-Landau theory for a single straight 
vortex line from the extrema of the energy density [2]. 
In an external magnetic field, there exist triangular and 
various other regular arrays of vortex lattices and vari- 
ous phase transitions. In the core of each vortex line, the 
pair field p goes to zero over a distance £. If we want 
to sum over grand-canonical ensemble of fluctuating vor- 
tex lines of any shape in the partition function (13), the 
space dependence of p causes complications. These can 
be avoided by an approximation, in which the system is 
placed on a simple-cubic lattice of spacing a = at;, with 
a of the order of unity, and a fixed value p — po given by 
Eq. (7). Thus we replace the partial partition function 
(17) approximately by 



Z 2 [po]= X! exp 

{1;V1=0} 



47r 2 p 2 1 a 



]Tl(xK oe (x-x')l(x') 



(18) 



The sum runs over the discrete versions of the vor- 
tex density in (12). These are integer- valued vectors 



l(x) = (Zi(x),Z 2 (x),Z 3 (x)) which satisfy V • l(x) = 0, 
where V denotes the lattice derivative. This condition 
restricts the sum over all l(x)-configurations in (18) to 
all non-selfbacktracking integer-valued closed loops. The 
function 
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«m(x) =]J / 



d 3 (aki 



i(kixi+k 2 x 2 +k 3 x 3 ) 



(2tt) 3 2J2 3 =1 (1- cos ah) +a 2 m 2 
- J d S e-^+ m2 >I Xl {2s)I X2 (2s)I X3 (2 S ). (19) 

is the lattice Yukawa potential [14]. 

The lattice partition function (18) is known to have a 
second-oder phase transition in the universality class of 
the XY-model. This can be seen by a comparison with 
the Villain approximation [15] to the XY model, whose 
partition function is a lattice version of 

z v [ P ] = Jve J ve v T exp - h - J d 3 x (ve-e v T y 

After integrating out 0(x), this becomes 

Z v [p] = Dct- 1/2 (-V 2 ) J V9 v T cxp (-1 J d 3 x8 v T 

(20) 

and we can replace B v2 by V x ^(-V 2 )" 1 ^ x 8 V T ) 2 = 
j v (— V 2 ) -1 ^. By taking this expression to a simple-cubic 
lattice we obtain the partition function (18), but with p^a 
replaced by (3y = ba, and the Yukawa potential Vp oe (x) 
replaced by the Coulomb potential i>o(x). 

The partition function (18) has the same transition at 
roughly the same place as its local approximation 



Z 2 [pa] ~ exp 

{1;V1=0} 



47r 2 p 2 a 



,(0)^l 2 (x) 



• (21) 



A similar approximation holds for the Villain model with 
v (x) instead of Vp oe (x) and p\a replaced by [3 V = ba. 

The Villain model is known to undergo a second-order 
phase transition of the VY-model type at (3y = r/3 with 
r w 1, where the vortex lines become infinitely long [17]. 
Thus we conclude that also the partition function (21) 
has a second-order phase transition of the VY-modcl 
type at p 2 Vp oe (0)a w v (0)/3. The potential (19) at the 
origin has the hopping expansion [16] 



v m (0) 



Hn 



V9 a ( a2m2 + 6 )" +1 

n=0,2,4 v ' 



H = 1,H 2 = 6, 



.(22) 



To lowest order, this yields the ratio v m (0)/vo(0) = 
l/(m 2 /6 + 1). A more accurate fit to the ratio 
v m {0)/vo(0) which is good up to m 2 w 10 (thus com- 
prising all interesting K-values since m 2 is of the order 
of 3/k 2 ) is l/(crm 2 /6 + 1) with a w 1.38. Hence the 
transition takes place at 



3 



Po a 



(aa 2 p 2 e 2 /6+l) 



or po 



V3a V 1 - orae 2 /18" 



(23) 



This transition will persist until p reaches the value 
Pi = 2c/3g of Eq. (9). From there on, the transition 
will no longer be of the XY-mode\ type but occur dis- 
continuously as a first-order transition. 

Replacing in (23) a by a£i of Eq. (10), and po by pi, we 
find the e quat ion for the mean-field Ginzburg parameter 
kmf = \/g/e 2 - 



I^MF "h ® ® 



0. 



(24) 



Inserting a « 1.38 and choosing a s=s r s=s f , the solution 
of this yields the tricritical value 



k>mf « 0.81/V2. 



(25) 



In spite of the roughness of the approximations, this re- 
sult is very close to the value 0.8/V2 derived from the 
dual theory in [11]. The approximation has three uncer- 
tainties. First, the identification of the effective lattice 
spacing a — a£ with a w 1; second the associated neglect 
of the x-dependence of p and its fluctuations, and third 
the localization of the critical point of the A^V-modcl 
type transition in Eq. (23) setting q 1. By modifying 
slightly the parameters q and r we can, of course, obtain 
complete agreement with the previous result.. 

Our goal has been achieved: We have shown the ex- 
istence of a tricritical point in a superconductor di- 
rectly within the fluctuating Ginzburg-Landau theory, 
by taking the vortex fluctuations into account. For 
k > 0.76/V2 these give rise to an XY- model type second- 
order transition before the cubic term becomes relevant. 
The cubic term causes a discontinuous transition only for 
k < 0.76/V2. 
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